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VARIATIONS IN THE SPATIAL DISTRIBUTION OF 11 MICRON 
RADIATION FROM OMICRON CETI* 

E.C. Sutton, J.W.V. Storey, and C.H. Townes 
Department of Physics 
University of California, Berkeley 



The spatial distribution of 11 micron radiation from o Ceti has 
been observed at various phases of its light cycle using a stellar 
interferometer. Changes have been seen which can be attributed to 
variation in the strength of thermal emission from circumsteTlar dust 
relative to the stellar continuum at 11 microns. These changes are 
shown to be cor:;elated with the changes ir, luminosity of o Ceti in 
such a way that dust grain emission at 11 microns was increased more 
than the continuum during the period of maximum luminosity. The de- 
gree of the change in dust grain emission implies that the maximum 
dust temperature is in the range of 500 K to 700 K during minimum 
-stellar luminosity. 


* This work was supported in part by the National Aeronautics and S’ jce 
Administration under grant NgL-DS^^0U3^^and in part by the National 
Science Foundation under grant AST 77-12256. 
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I INTRODUCTION 

Long-period variable stars cownonly are surrounded by circumstellar 
dust shells. Variations in, the luminosities of these stars will greatly 
change the radiative environment of the dust grains; thus it is expected 
that variations will be seen in the thermal emission from these grains. 
Conventional techniques which might detect such variations include mea- 
surements of broad-band infrared colors and also narrow-band spectro- 
photometry in the 8-14p region. Forrest, Gillett, and Stein (1975) have 
reported such measurements on o Ceti and several other long-period 
variables. They found no evidence for large changes in the strength of 
’ emissiW f r^tiS T. The work reported in this letter 

shows that there are changes in the dust grain emission which can be 
seen using rather different techniques. These variations are consistent 
with the limits set by Forrest, Gillett, and Stein. 

Interferometric techniques in the infrared may be used to determine 
directly the r-^lative brightness of circumstellar emission with respect 
to the stellar continuum and also the sizes of these regions of circum- 
stellar emission (Sutton et al. 1977; McCarthy, Low, and Howell, 1977). 
For this. Wavelengths near 11 microns are of particular importance due 
to the spectral feature attributed to silicate grains. Spatial resolu- 
tions of existing instruments at this wavelength are of the order of 1" 
which is an appropr^iate scale for examining many circumstellar dust 
shells. It is thus possible to examine variations in thermal emission 
from dust around long-period variables in more detail than previous 
techniques have allowed. 

II OBSERVATIONS 

The interferometer used in these observations makes use of the 
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McMath auxiliary telescopes at Kitt Peak 'National Observatory.** The 
equipment and techniques used have been described by Sutton et aK (1977) 
During these observations the sensitivities of the heterodyne receivers 
were approximately 3 x 10“ watts/zllz' for an infrared bandwidth of 
0.1 cm'^ at 11 . 106p. 

Observations were made of o Ceti from February 1977 through February 
1978, which represents slightly more than one cycle of the visible 
light curve. No observations Were possible for several months before 
and after May 1977 during which time o Ceti was in the daytime sky. 

The overall instrumental cal ibration, original ly based on labora- 
tory measurements, was verified by observing the star a Her which has 
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1976). Thus tt Her should behave as a point source since its diameter 
is at least an order of magnitude smaller than 0.4", which was the 
finest lobe spacing used. Three observations were made yielding an 
average fringe visibility of 0.95 ± 0.08. Since this result was con- 
sistent with unit visibility, no further correction was applied to the 
visibility scale. As a secondary standard, a Ori was observed with 
maximum resolution within a few days of each observation of o Ceti. 

The observations of a Ori were internally consistent. 

Figure 1 presents. the measured fringe visibilities for o Ceti. 

The phase quoted is with respect to the visible light curve with = 0 
representing maximum light. Note that Figure Ic contains data from two 
separate cycles but at nearly the same phase. Figure H shows a compar- 
ison set of observations on a Ori which indicate the flatness of the 
instrumental response. 

** Operated by the Association of Universities for Research in Astronoii\y, 
Tnc. , under contract with the National Science Foundation. 
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The observations have been fitted with a two-component model of 
the brightness distribution in which one component is spatially unre- 
solved and emits a fraction, f, of the flux at lip. The other compo- 
nent, which contains the remaining 1-f of the flux, was assumed to have 
a Gaussian spatial distribution with a full-width to 1/e intensity 
given by the parameter 6g. The curves drawn in Figure 1 are the best 
fits to the data using this model, and the values of f and 6^ used to 
make these fits are listed in Table 1. 

■ ' I 

This two-component model was chosen because it was known that 
there would be contributions to the lip flux both from the stellar 
eontinuum- and" from dual; graiTr emi 551 07 \t The “dfameter" of' i;he phofo-^^" " 
sphere of o Ceti is sufficiently small that the region of stellar con- 
tinuum emission should appear point-like at these resolutions (Labeyrie 
et al., 1977). This is demonstrated by the flatness of the visibility 
curves in Figure 1 at high resolution. The component representing dust 
grain emission was arbitrarily assumed to have a simple Gaussian spatial 
distribution. None of the results presented here are greatly affected 
by thr detailed shape of this distribution. 

The most noticeable variation shown in Table 1 is in the parameter 
f. Since 1-f is the fraction of the 11 micron flux contributed by dust 
emission, it is seen that this dust emission is strongest relative to 
the stellar continuum near a phase of 0.10 or slightly later. This 
phase also corresponds to maximum stellar luminosity indicating that 
the dust grain temperatures are raised by the increased stellar radiation 

The data in Table 1 also suggest some variation in the spatial ex- 
tent of the region of cirGUjmstellar dust emission. The diameter, o , 
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seems to increase together with the strength of dust emission, indicating 
that it also may be related to the increased stellar- luminosity. 

Changes in stellar luminosity may be caused by changes in either 
effective temperature or size. The effective temperature of o Ceti is 
thought to undergo a periodic variation from about 2300 K to 2800 K 
(Strecker, 1973). This alone predicts a factor of 2.2 change in the 
total stellar luminosity. Small diameter changes may also be present, 
but Strecker calculates that the luminosity change is still about a 
factor of 2.2. Although this calculation is subject to a number of 
uncertainties, it is in good agreement with the factor of 2.5 variation 
in bolometric magnitude observed by Pettit and Nicholson (1933). The 
mWii^ifli weTT ^as the maxima of the 3 . 5u and 1 1 n 1 d Qht . .. - 

curves fall at a phase of approximately 0.15 after visible maximum. 

Thus the changes in luminosity are of the correct phase to explain the 
variations in dust grain emission. 

The strength of dust emission at 11 microns relative to the 
stellar continuum is seen from Table 1 to vary by a factor of 1.8. As 
the stellar temperature varies from 2300 K to 2800 K, the 11 micron 
stellar continuum will itself increase by a factor of 1.3. Thus, the 
overall change in dust emission is a factor of 2.3. This is to be 
compared with the calculated effect of changes in stellar luminosity 
on the strength of dust emission. 

The expected magnitude of the variation in dust emission depends 
on the temperature of the grains. It is likely that emission is seen 
from dust at a variety of distances from the star and hence at a va- 
riety of temperatures. The maximum temperature present, is limited hy 
the condensation temperature of the grains, which is ap()roximately 
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1000 K for silicate materials (Sal peter, 1977). The actual temperature 
of the hottest grains may be considerably lower if grain formation is 
inhibited due to difficulties in nucleation. The lowest temperature of 
interest is 250 K since dust grains cooler than this will not radiate 
strongly at llu. Grains at 1000 K and 250 K v,‘il1 change in brightness 
by factors of 1.4 and 2.6 respectively when subjected to a factor of 
2.2 change in stellar luminosity. The expected variation in brightness 
for the entire dust shell should lie between these extremes and thus 
is consistent with the observed factor of 2.3 variation in dust emission. 
The observations further indicate that the radiation must come predom- 
inantly from grains near the cooler end of this range of temperatures. 

Detailed calculations have been made for the brightness variation 

of an optically thin dust region using several models. It was assumed 

that the dust region has an inner radius characterized by some maximum 

temperature for grains in radiative equil ibrium and that the dust 

density follows a power-law distribution, r"*^, outside of this radius. 

For grains with a l/\ emissivity spectrum and a maximum temp.erature of 

1000 K a power-law given by n = 1.5 is needed to produce the observed 

variation. This particular result does not depend strongly on the dust 

spectrum. The more likely density dependence n = 2 requires a Tower 

maximum temperature between 500 K and 700 K. Models with grain density 

-15 

dropping off slower than r ' or with still lower temperature dust are 
ruled out in all cases since the emphasis on cooler grains would produce 
variations in the llu shell brightness with stellar phase much larger 
than those observed. 

The observed shell diameter allows a rough deteniiiriation to Ikj ttujf)(‘ 
of the grain emissivity spectrum. The temperature of grains at the ififier 
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boundary of the shell has been shovm to be between 500 K and 700 K. 

The grain einissivity spectrum determines the distance from the star at 
which the grains will be in equilibrium at this temperature. The ob- 
served diameter of 0.7" fits the brightness distribution of a model 
with an emissivity spectrum of approximately 1/X. 

The diameter of the shell appears to change somewhat with stellar 
phase. Such a change is the result of the fact that the cool grains 
change lly brightness faster than hotter grains as the stellar . luminos- 
ity increases. Thus emission at larger distances from the star is em- 
phasized a^ the dust is heated, The model with a maximum grain tempera- 

_2 _i" 

ture of 500 K to 700 K, an r density dependence, and a X grain 

emissivity spectrum predicts a 20% increase in diameter from minimum 
to maximum stellar luminosity. This increase is consistent with the 
observed variation within the rather large errors. This apparent size 
change need not involve any physical redistribution of the dust itself. ' 
Such redistribution is possible although unlikely within the timescales 
involved here. 

The optical deeth of the dust shell can be derived from the bright- 
ness of the circumstel lar shell, but it is dependent on all of the param- 
eters of the model. The result for the model indicated here is an 
optical depth at 11 microns of approximately 0.02. 

Since the properties of this dust model depend rather strongly on 
the observed variation in dust brightness, it is important to consider 
other more conventional measurements which may be able to detect this 
variation. For example, the 8-14ii spectrum of o Ceti should show a 
more pronounced silicate emission feature at phases near maximum lumi- 
nosity. Such a variation is suggested in the data of Forrest, Gillett, 
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and Stein (1975). However, this should be interpreted with caution 
since their single spectrum near maximum light predated their other 
spectra by several years. In addition, it is possible that such spec- 
tra may have been distorted due to phase-dependent changes in the 
strength of SiO absorption bands at 8p such as those changes detected 
by Hinkle et al. (1976) in the first overtone bands of SiO at 4ii. If 
so, as maximum light was approached the 8g absorption may have weakened 
at the same time as silicate emission near Hu was getting stronger. 

Thus the shape of the spectrum over this narrow region may have remained 
relatively unchanged. Nevertheless, further investigation of the 8-14u 
spectra of long-period variables at various phases may prove fruitful. 

Broadband infrared color measurements provide additibhaTnhfb^^^ 
tion. The observations of Forrest, Gillett, and Stein show a variation 
of approximately in CHu3 - C8.4u3 with the Hu band being re- 

latively brightest at phases between 0.0 and 0.Z5. As with the spectro- 
photometric measurements the effect of SiO absorption near Bp may be 
such as to mask an'inherently much larger variation in tHu3 - C8.4u3. 
Similarly the quantity Clip) - [3,5p] is contaminated by the H^O absorp- 
tion band at 2.7p. Forrest, Gillett, and Stein observed virtually no 
variation in this color as a function of phase. Yet considering the 
HS*” amplitude of the 2.7u light curve observed by Maran et al . (1977) 
and the width of the absorption in this region as shown by the data of 
Merrill and Stein (1976), it is to be expected that there would be a 
variation in [3.5p3 of approximately 0.3^ solely due to changes in H 2 O 
absorption. This effect masks variations in CHiH - fS-Siil which may 
be present due to increased dust emission. Thus the broadband photo- 
metric measurements are consistent with the observed increased dust 


emission at phases near 0,10 although by themselves they are not 
sufficient to demonstrate such an increase. 

IV SUMMARY 

The observed variations in the 11 micron brightness distribution 
of 0 Ceti indicate an increase in both the relative brightness and 
the apparent size of the circumstellar dust shell occurring during 
the period of maximum stellar luminosity. These variations can be 
understood in terms of changes in the temperature of the dust grains 
produced by the known variation in luminosity of o Ceti. 



The authors especially wish to thank L. Greenberg and A. Harris for 
their assistance in making these observations and A. Betz for his help 
with the equipment. 
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TABLE 1 

PARAMETERS iJSED TO FIT FRINGE VISIBILITY MEASUREMENTS 


P'otire 

Phase 


f 





la 

0.75 

0.54 

± 0.02 

*■ 


lb 

0.90 

0.54 

O 

o 

ro 

o'“!62 ± 

0.08 

Ic 

0.11 

0.40 

+ 0.01 

o''!77 ± 

0.14 

1-/4 . . . 


Q-AR 

+ 0.02 

fi^7Q L 

on? 


* Undetermined due to the lack of measurenents at low spatial resolution. 


FIGURE CAPTIONS 


Figure 1: Fringe visibility measurements of o Ceti at 4 different 

phases of its Tight cycle. Curves drawn are fits using 
the parameters of Table 1 as discussed in the text. The 
graphs are in order of increasing phase centered about 
^ = 0 . , 

Figure "2 r^^^ompaH^rr firings visibiTity -me3suremen^fes. fQE a Qrio 
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